Pressure-induced s-band ferromagnetism in alkali metals 
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First-principles density-functional-theory calculations show that compression of alkali metals sta- 
bilizes open structures with localized interstitial electrons which may exhibit a Stoner-type instability 
towards ferromagnetism. We find ferromagnetic phases of the lithium-IV-type, simple cubic, and 
simple hexagonal structures in the heavier alkali metals, which may be described as s-band ferro- 
magnets. We predict that the most stable phases of potassium at low temperatures and pressures 
around 20 GPa are ferromagnets. 

PACS numbers: 71.15.Nc,75.10.Lp,62.50.-p,61.66.-f 



o 



> 

rn 
O 



X 



Among the elements, bulk ferromagnetism is found 
only in the first row of the transition metals and in the 
lanthanides, where partially filled 3d and 4/ electronic 
shells of strongly-localized orbitals are present. Magnetic 
order is normally reduced and eventually destroyed by 
the application of pressure because it tends to delocalize 
electronic states. The alkali metals are archetypal nearly- 
free-electron materials in which each atom contributes a 
single valence electron to a gas whose effective interac- 
tion with the ionic cores is weak. They might there- 
fore be thought of as the least likely elements in which 
to find bulk ferromagnetism. The alkali metals, lithium 
(Li), sodium (Na), potassium (K), rubidium (Rb), cesium 
(Cs), and probably francium (Fr) (the most unstable 
naturally occurring element) adopt body-centered-cubic 
(bcc) phases under ambient conditions which, on com- 
pression, transform to face-centered-cubic (fee) phases 
[l|. Diamond-anvil-cell experiments have shown, how- 
ever, that they adopt more open structures at higher 
pressures Density- functional theory (DFT) calcula- 
tions have reproduced the stability of the experimentally 
observed phases and have given insights into their elec- 
tronic structures • The most surprising result is that 
a large amount of valence charge in the open structures 
resides within the interstitial regions rather than in close 
proximity to the ions. This corresponds to the formation 
of "electrides" in which the interstitial electrons form the 
anions d-Q. The open structures correspond to well- 
packed ionic solids when both the alkali metal ions and 
the centers of the interstitial electronic charges are des- 
ignated as ionic positions 

The band structures of alkali metals deviate substan- 
tially from nearly-free-electron behavior under applied 
pressure The occupied valence bands become 

flatter than the corresponding nearly-free-electron ones, 
which narrows the occupied valence bandwidth. Ashcroft 
and coworkers 0, 0, 0] have attributed this phenomenon 



to the interaction of the valence electrons with the rela- 
tively incompressible ionic cores that occupy an increas- 
ingly large fraction of the total volume as pressure is 
increased. Under ambient conditions the effective in- 
teraction between the valence electrons and ionic cores 
is weak, but it becomes strongly repulsive under pres- 
sure and forces valence electrons to occupy interstitial 
positions 0, 01 ■ The interstitial regions in close-packed 
structures are numerous, but small, and "cutting up" the 
valence charge into small regions increases the kinetic en- 
ergy. The kinetic energy can, however, be reduced by 
adopting more open structures which have less numerous 
but larger interstitial regions in which to accommodate 
the valence electrons. This effect evidently overcomes the 
concomitant increase in core-core repulsion energy. 

The Fermi surface/Bragg plane mechanism for stabilis- 
ing structures [lo| involves the action of the Fourier com- 
ponents of the lattice potential on the degenerate electron 
orbitals on the Bragg planes in reciprocal space. The to- 
tal energy is lowered if a structure is adopted in which 
Bragg planes graze the Fermi surface, because the occu- 
pied orbitals just inside the Bragg plane are lowered in 
energy while those of the unoccupied orbitals just out- 
side are raised. This mechanism can lead to structural 
instabilities and it is likely to be involved in determining 
the details of the high-pressure structures of many of the 
alkalis 11, 12[. The Fermi surface/Bragg plane mecha- 



nism normally operates when g{E) is large in the region 
around the Fermi energy Ep, but an alternative insta- 
bility comes into play when g{E) is also strongly peaked 
around E-p. In this case the energy may be lowered by 
transferring electrons from one spin channel to the other, 
resulting in ferromagnetism (FM). This is the Stoner in- 
stability [13[. The key quantities in Stoner's theory are 
g{Ep) and the effective interaction / between the up and 
down spin densities in the unpolarized state, with the 
instability to FM occurring when g{Ep)I > 1. 
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Motivated by the above analysis, we have investigated 
magnetic ordering in the alkahs under high-pressures. We 
searched for low-enthalpy structures using ab initio ran- 
dom structure searching (AIRSS) which has been 
successfully applied to systems as diverse as metals un- 
der high pressures 0, Q and molecular solids flil 17 [. 
We used the CASTEP plane-wave DFT code and 
the Perdew-Burke-Ernzerhof (PBE) Generalized Gradi- 
ent Approximation (GGA) density functional [3], and 
we also present results obtained with a Thomas-Fermi 
screened exchange functional (SX) (2^] and the local spin 
density approximation (LSDA). We used ultrasoft pseu- 
dopotentials 2lj, treating all three electrons explicitly 
for Li, and nine electrons for the other alkalis. Brillouin 
zone integration grids of spacing 27rx0.05 were used 
for the searches and the low-enthalpy structures were fur- 
ther relaxed using a finer grid spacing of 27rx 0.03 A-i 
which gave very accurate enthalpy differences between 
the phases. Large plane wave basis set energy cutoffs 
were used. The phonon calculations were performed us- 
ing a finite-displacement method and 64-atom supercells. 

We performed spin-polarized calculations, starting 
some relaxations in a high spin state with an average spin 
density of one electron per atom and others with zero av- 
erage spin density. In each calculation the spin density 
was allowed to evolve freely as the structure was relaxed. 
We also performed calculations without spin polariza- 
tion. The use of a wide variety of starting structures and 
spin states was found to be important in allowing many 
different spin and atomic configurations to be accessed. 
The main searches were performed at pressures around 
those at which experiments show that the fee structures 
transform to more open phases. We performed searches 
with unit cells containing up to eight atoms, relaxing a 
total of about 1500 structures. Calculations for other 
known structures of the alkalis which have more than 
eight atoms per cell were also performed, namely a pe- 
riodic structure of space group 1 4/ mem with 56 atoms 
which is a good analogue of the K-III (Rb-IV) incom- 
mensurate host-guest structure, and the 84-atom C222i 
structure of Cs-III and the related 52-atom structure of 
Rb-III [1]. 

We find good agreement with the experimentally ob- 
served phase transitions. Our calculated coexistence 
pressure for K-fcc and K-IA/mcm of 19.8 GPa is close 
to the experimental transition pressure of 23 GPa T]. 
We note that Marques et al. [23j have reported finding 
the K-/4/rncm phase in some high- pressure experiments 
and the K-hP4 phase in others, which is consistent with 
our finding that K-hP4 is only slightly less stable than 
K-I4/mcm. Our coexistence pressure for the Rb-fcc and 
Rb-C222i phases of 14.2 GPa is in excellent agreement 
with the experimental transition pressure of 14 GPa [ij . 
We obtain a coexistence pressure for Cs-fcc and Cs-C222i 
of 4.8 GPa, compared with the experimental transition 
pressure of 4.2 GPa 1]. We find Cs-IV to be stable from 



4.8-10 GPa, in excellent agreement with experiment [l[. 
We find no region of stability for the observed Gs-III 
structure which in our calculations only becomes more 
stable than the fee phase at 4.9 GPa, although this is a 
small discrepancy. In general we find that our calculated 
coexistence pressures can be brought into agreement with 
the experimental data by rigidly shifting the enthalpy 
curves by less than 5 meV, which indicates the high level 
of accuracy of our calculations. 

We found M3d, simple cubic (sc) and simple hexagonal 
(sh) phases of K, Rb, Cs and Fr with strong FM ordering, 
and a weakly FM Cs-fcc phase, see Fig. [TJ The enthalpy 
reductions due to the formation of FM moments in the 
sc phases are similar to those in the corresponding sh 
phases. The sh phases are more stable than the sc phases 
at lower pressures in K and Rb, and more stable at all 
pressures in Cs. We did not find any spin-polarized states 
of Li or Na in our fully- converged calculations, but we can 
obtain them by reducing the number of k-points. This 
shows that Li and Na are close to a FM instability and 
demonstrates the importance of carefully studying the 
convergence with respect to the k-point sampling, as we 
have done. We predict that the FM K-M3d phase to be 
the most stable in the range 18.5-20 GPa and the FM 
K-sc phase to be the most stable in the range 20-22 GPa. 
The FM ordering leads to an energy gain of a few tens 
of meV per atom. When the centers of the interstitial 
electronic charges are designated as ionic positions, sc 
becomes the CsCl structure, sh the MgB2 structure, and 
/43d the Th3P4 structure with the Cs, Mg and P 
sites being those of the cations and the CI, B and Th 
sites being those of the interstitial electrons. 

To explore the sensitivity to the density functional we 
also performed calculations using a Thomas-Fermi SX 
functional j20| and the LSDA. For the SX calculations we 
used a screening wave vector of kg — 0.764 a.u., which 
corresponds to a Wigner-Seitz radius of rg = 4.2 a.u. 
This is larger than the value of = 3.4 a.u. obtained 
from the average valence charge density of FM K-sc at 
20 GPa, but the results are rather insensitive to reason- 
able variations in ks as it is proportional to the one-sixth 
power of the average charge density. Using the SX func- 
tional instead of PBE-GGA stabilizes FM K-sc over K- 
fcc by about 15 meV per atom. We expect that the sc 
phases of each of the alkalis would be further stabilized 
with respect to the fee phases when calculated with the 
SX functional. We did not find magnetic ordering in 
K when using the LSDA functional which predicts the 
K-sc phase to be stable above about 20 GPa, in disagree- 
ment with experiment. Given that the pressures of the 
observed transitions calculated within PBE-GGA are in 
very good agreement with experiment, we are inclined to 
believe that it gives a better description of compressed 
alkalis than the LSDA or SX functionals. 

Calculations of the harmonic vibrational modes of the 
K-fcc, K-sc and FM K-sc phases showed them to be dy- 
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namically stable. The phonon modes of the sc phases are 
substantially softer than those of K-fcc, and including 
the zero-point enthalpy stabilizes the K-sc and FM K-sc 
phases over K-fcc by about 25 meV per atom at 20 GPa. 
The stability of K-sc and FM K-sc over K-fcc is slightly 
increased by including the vibrational contribution to the 
free energy at a temperature of 300 K, but the magnetic 
moment will be reduced by disordering of the spins. 

The valence electron density of states of the paramag- 
netic and FM states of K-sc at 20 GPa are shown in Fig. 
[2j In the paramagnetic system, Ep is close to the top of 
a large peak in g{E) and the system is ripe for a Stoner 
instability. In the FM K-sc phase Ep falls just above the 
peak for the majority-spin band and well below the peak 
for the minority spin band. The FM K-sc state is about 
35 meV per atom more stable than the paramagnetic K- 
sc phase, which is sufHcient to make FM K-sc the most 
stable phase in the pressure range 20-22 GPa. The value 
of g{Ep) for non-spin-polarized K-sc at 20 GPa is 1.4 per 
eV per atom, so that the Stoner criterion is satisfied for 
/ > 0.7 eV, which is similar to the values deduced for 
FM in transition metals j24| . 

The spin densities of FM K-sc and FM K-sh at 20 
GPa are shown in Fig. [31 A large blob of spin-polarized 
charge density resides on the interstitial site at the cen- 
ter of the cube of K-sc, while the spin polarization on the 
atomic sites is small. The spin polarization on the atoms 
is also small in FM K-sh, and the interstitial spin den- 
sity is more diffuse. The FM K-sc structure achieves its 
maximum spin moment of about 0.72 electrons per atom 
at 22 GPa, while the maximum spin moment in FM K- 
sh of about 0.62 electrons per atom occurs at 20 GPa, 
and the maximum spin moment of K-I43d is about 0.4 
electrons per atom. The FM phases might be described 
by a Hubbard-like model [1^ using tight-binding s or- 
bitals centered on the interstitial regions. The magnetic 
state would then be described as s-band FM, which is 
very different from the d oi f band FM observed in the 
transition metal and lanthanide elements. 

We are not aware of any experimental evidence for 
magnetic ordering in bulk alkali metals, although there 
is evidence of it in low-dimensional systems. Theoretical 
work by Overhauser suggested that charge-density 
and spin-density-wave instabilities might occur in alkali 
metals at low pressures, but these ideas have not been 
widely accepted. A DFT study by Zabala et al. 27 1 
found spontaneous magnetization in "Na wires" mod- 
elled by jellium, and Bergara et al. [2^ found instabilities 
to FM in atomically-thin Li and Na wires. FM ground 
states are, however, forbidden in strictly one-dimensional 
systems by the Lieb-Mattis theorem [29| and, in quasi- 
onc-dimensional systems, quantum fluctuations tend to 
suppress FM. Experimental observations of FM in K clus- 
ters incorporated within a zeolite [s^l and anti-FM in K 



clusters in a nanographite-based host 3l| have also been 
reported. We have used the AIRSS method to search for 
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FIG. 1: (color online). Enthalpy-pressure curves for the FM 
and paramagnetic phases of K, Rb, and Cs. Differences in 
enthalpy from the fee phase are plotted. FM phases are shown 
as solid lines and paramagnetic phases as dashed lines. 
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FIG. 2: (color online). The electronic density of states, g{E), 
of K-sc at 20 GPa in states per eV per atom. Data for the 
paramagnetic phase is shown in brown and for FM K-sc in 
grey. The Fermi energy is shown as a vertical dotted line. 



structures of small unsupported K clusters, finding that 
some weakly anti-FM states are energetically favorable 
and that quite strongly FM phases with moments up to 
about 0.5 electrons per atom occur at higher energies. 
Electride formation in alkali metal clusters at ambient 
pressure is highly unlikely because the bonds are too long, 
and the magnetic ordering must produced by some other 
mechanism. It might be possible to further stabilize FM 
phases of bulk alkali metals by alloying them with other 
alkalis or other species. Low-temperature experiments 
looking for magnetism in compressed alkalis are needed 
to test our predictions. 

In summary, we predict mechanically-stable FM 
phases of the heavier alkali metals to have low enthalpies 
at pressures just above the stability range of the fee 
phases. FM phases with the Li-IV (lASd), sc, and sh 
structures can be described as s-band electride ferromag- 
nets. The FM K-I43d and K-sc phases are predicted to 
be the most stable at low temperatures and pressures 
around 20 GPa. 

The authors were supported by the Engineering and 
Physical Sciences Research Council (EPSRC) of the UK. 




FIG. 3: (color online). The spin density of FM K-sc (top) and 
FM K-sh (bottom) at 20 GPa. The purple spheres represent 
the atoms and the spin density is shown in green. 
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